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Abstract 

The performance of the active air solar heating system 
installed in Module 6 of the Alberta Home Heating Research 
Facility was evaluated for the September 1981 to March 1982 
heating season. The operating characteristics and energy 
flows for the system were examined and the performance 
parameters calculated. 

The solar contribution assessed by direct and indirect 
methods ranged from 11.3 to 17.8 %. The indirect methods 
yielded values that were close to the measured solar 
eontrabution. 

The solar collector was constructed from six commercial 
collectors and its performance was calculated on an 
instantaneous, hourly, daily and monthly basis. For this 
solaricollectorp;ethéaFr(ra)e ranged from 0.53 to 0.56, and 
all values were higher than the performance test value of 
0.5. The collector heat loss coefficient FrUt ranged from 
(ASLO 2.461 W/m? =2e, The monthly value Fru of 2.61 w/m?—-C 
is closer to the manufacturer's performance test value of 
2.9 W/m?-°C and is probably more representative of the 
collector heat losses. 

Energy flows for the components of the system were 
evaluated. The 25m rock bed storage bin was fully charged 
during periods in September 1981 when heating requirements 
were lows About 23.6 4vot the 3525.3 Md of mrotal venergy 


stored was transferred from the bin. 
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Nomenclature 


Notation 

Ac Transparent frontal area for a flat-plate collector, 
[m? ] 

Ag Gross collector area, [m?] 

COL Operation time for heat-collection mode, [min] 

ep Specificyhéateof the’ transfer fluid, air, [J/kgq-°C] 

EH Operation time for electric heating mode, [min] 


DTCcol Hourly sum of temperature difference across 
collectors in heat-collection mode, [°C-s] 

DTChfc Hourly sum of temperature difference across 
collectors in heat-from-collector mode, [°C-s] 

DTScol Hourly sum of temperature difference across storage 
bin during heat-collection mode, [°C-s] 

DTSeh Hourly sum temperature difference across storage 
bin during electrical heating mode, [°C-s] 

DTShfs Hourly sum of temperature difference across storage 
bin in heat-from-storage mode, [°C-s] 

are Instantaneous temperature difference measured 
Heross coéolbéctors pf ITiesTiige[ Sc) 

agTs Instantaneous temperature difference measured 
Mecossustoragqe sbinverst-Tsbyel cd 

£ Monthly solar contrcrbutionvassessed*by Gk =chart 


prediction, dimensionless 


Fua. Power correction factor, dimensionless 
Fr Solar collector heat removal factor, dimensionless 
HFC Operation time for heat-from-collector mode, [min] 
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HFS Operation time for heat-from-storage mode, [min] 

Ht Total solar radiation on the tilted surface, [W/m] 

I Total no. of hours corresponding to each operation 
mode in a day, [Hr] 

K Total no. of hours corresponding to each operation 
mode in a month, [Hr] 

m Mass flow rate of air, [kg/s] 

Qas Actual monthly solar energy supplied from the 
thermal system to space heating, [J] 

Qbhl Hourly energy losses from the module through 
conduction, [WwW] 

Oecol Hourly solar energy gained from collector during 
heat-collection mode, [J] 

QOchfc Hourly solar energy gained from collector during 


heat-from-collector mode, [J] 


Qeh Hourly energy supplied during electric heating 
mode, [J] 

Qeh5 Hourly electric power usage in module 5, [W] 

Qeh6 Hourly electric power usage in module 6, [W] 

Qht Hourly total energy supplied or required for 


space heating, [J] 
Qmeh Monthly electric energy usage, [J] 
Qmsun Monthly solar energy gained from collector, [J] 
Qmas Monthly actual solar energy supplied from the 
thermal system to space heating, [J] 
Qmsa Monthly total active solar energy required for 


space heating, [W] 
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Qscol Hourly energy entering the storage bin from 
eolvector,e {a ] 

Qseh Hourly energy supplied from storage bin for space 
heating during electric heating mode, [J] 

Oshfs Hourly energy supplied from storage bin for space 
heating during heat-from-storage mode, [J] 

Qsloss Hourly energy losses from storage bin by 
air leakage or conduction through bin walls, [J] 

Qsto Hourly energy stored or accumlated in the storage 
ba ee fo] 

Qsysl Monthly energy transfers or losses from the 
thermal system, [J] 

Ta Ambient air temperature, [°C] 

aie Temperature of transfer fluid entering the 


collector? lpecd 


Tfe Temperature of transfer fluid leaving the collector, 
neal 
iu) Average temperature of the collector absorbing 


surtace , Miecs 
Tr Room temperature, [°C] 


Tsb Temperature of transfer fluid at the bottom of 


storage bin, [°C] 


Fst Temperature of transfer fluid at the top of storage 
banger [fed 

UA Building heat loss coefficient (UA Factor), [W/°C] 

Ed Daily collector efficiency, dimensionless 

Eh Hourly collector efficiency, dimensionless 
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£1 Steady state instantaneous collector efficiency, 


dimensionless 
—m Monthly collector efficiency, dimensionless 
Ls Monthly fraction of energy transfers or 


losses from the storage bin, dimensionless 

Lsd Monthly fraction of energy transfers or 
losses from the storage bin and the duct section, 
dimensionless 

Lsys Monthly fraction of energy transfers or 
losses from the thermal system, dimensionless 

nbhl Monthly solar contribution assessed by building 
heat loss method, dimensionless 

neh Monthly solar contribution assessed by power 
records method, dimensionless 

Nr Monthly solar contribution assessed by solar 
required method, dimensionless 

nS Monthly solar contribution assessed by solar 
supply method, dimensionless 

Py Monthly power ratio for module 5 and 6, 
dimensionless 

Pya Monthly actual power ratio for module 5 and 6, 
dimensionless 

a Absorptance of the collector absorber surface for 
solar radiation, dimensionless 

T Transmittance. of the solar collector covemmplate, 


dimensionless 
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Effective transmittance-absorptance product at 
normal incidence, dimensionless 

Effective transmittance-absorptance product, 
dimensionless 

Solar collector heat transfer loss coefficient, 
Lw/m?=2C4 


Air Density, [kg/m?] 
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1. Introduction 

The production rate of petroleum liquid for the North 
American continent appears to have reached a maximum and may 
now be decreasing. Foreign sources of petroleum liquid are 
the major supplementary sources of energy and these also may 
be exhausted in the foreseeable future. In order to ensure a 
continuous supply of energy, the adoption of conservation 
measures and the developing and deploying other energy 
sources are necessary. AS a consequence, a number of 
researchers and practicing engineers changed their focus 
from non-renewable fossil fuel technologies to renewable 
power sources or technologies for reducing the consumption 
of fossil fuels. Many people in this field agree that solar 
energy is one of the most Significant of the renewable 
energy sources. The potential of solar energy assessed by 
the U.S. National Science Foundation in the testimony before 
the U.S. Senate Interior Committee in 1972 was that solar 
energy is the most promising of the unconventional energy 
sources potentially capable of meeting a significant portion 
of the future energy needs[28]. 

In recent years, residential energy consumption has 
become a very topical subject studied by many researchers. 
All of them tried to reduce auxiliary space heating 
requirements by a number of solar heating schemes, both 
active and passive. The result was many demonstration houses 
that modified the present house design in Canada and the 


United States. 
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The Alberta Home Heating Research Facility, which was 
constructed in the summer of 1979, is an experimential 
hesearchutaci lirtyreomsisting ofdsixréal neh. 3nxconQiemtsprigle 
story modules with full basements. These modules are located 
on the University of Alberta Farm ina Single East-West row 
and are designed to test domestic heating strategies ina 
northern climate. The modules are unoccupied. Two of the 
modules, module 5 and 6 were constructed firstly to compare 


the performance of two identical constructions under in situ 


conditions, and secondly to assess the operating 
characteristics of two active solar systems with different 
working fluids operating under the same meteorological 
conditions. However, because of budgetary constraints, only 
the active air heating system was installed. Details of this 
facility are presented in Chapter 2. 

Since the performance of a solar air heating system is 
a complex function of weather, system design, and collector 
area in relation to the load, a meaningful rating system 
also becomes very complex. A number of researchers such as 
Lunde[8],[9]; Klein, Beckman, Duffie[11]; Lof, Shaw, 
Oank[17]; Hill[23]; Gupta,Garg[14]; Hollands, Chineck, 
Chanarashekar[42]; and Close[6] proposed many methods of 
studying the monthly and annual performance of solar air 
heating systems. All of them emphasised the thermal 
performance of the collector and some considered both the 


performance of collector and storage bin. 
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Also, interactive computer programs for solar heating 
and cooling design have been developed by many designers. 
Two computer programs are in widespread use in connection 
with modeling or sizing of solar air heating system. The 
program called TRNSYS, which is for simulating the hour by 
hour dynamic thermal behavior of active solar heating and 
cooling systems, was developed by Klein, et al.[31]. Another 
program called F-chart, which is for simulating the 
long-term performance of solar heating or cooling system, 
was developed by Beckman, Duffie and Klein[30]. 

Although several methods for rating of solar heating 
systems have been proposed, it seems that many of them do 
not fully meet the requirements of technical and 
non-technical users of solar heating system. For the 
technical users, the solar collector performance, the 
storage bin performance and the solar space heating 
contribution are the main parameters to assess the 
performance of a solar heating system. The collector 
performance represents the heat delivery capability of the 
solar collector and is usually obtained from a graph of 
collector! efficiencyyplottedvasra functuoneetethen difference 
in inlet and ambient temperature divided by solar radiation 
onuthe. collectorméiditeds surfaces (Tfi-fad/Hte Thisheoldector 
efficiency curve can be applied to computer programs to 
assess the solar contribution. For the non-technical user, 
the solariheatung GContrabutdonai sothehmaineparameteérrito 


assess the thermal performance of a solar heating system. 
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Several national standards for rating solar collector 
have been developed, including those of the Air Conditioning 
and Refrigeration Institute[3], (ARI), the Solar Energy 
Industries Association[4], (SEIA), and the Solar Energy 
Research and Education Foundiation [5], (SEREF). All of 
these rating methods are based on the same collector testing 
method. This method, which was developed by the National 
Bureau of Standards and further refined by the American 
Society of Heating, Refrigerating, and Air Conditioning 
Engineers is now known as the ASHRAE 93-77 collector testing 
method[1]. However, the testing is done under cloudless sky 
conditions with uniform solar radiation and may not be 
representative of the dynamic behaviour of solar collectors 
in the real environment. 

The present study is mainly concerned with the 
performance of the active air heating system during the 
September 1981 to March 1982 heating season. The operating 
characteristics and energy flows for the system are examined 
and the performance parameters are calculated. Methods of 
assessing the solar contribution for the technical and 
non-technical solar heating users are presented in Chapter 
3. The collector performance examined on an instantaneous, 
hourly, daily, and monthly basis is presented in Chapter 4. 
The storage bin characteristics are also included in this 
chapter. 

In addition, the assessment of the system performance 


should include the effect of energy losses. Losses from the 


pitnedsibaed. ih ‘ite % sheds iiibasai | 
VOIPAR sede) oi wfiniad (edad taadie-e 

| teise@ webs? se bas, < (aga) ppeaperiaty: 
‘ap C0 . (sage) itaye tofsatBoud notseoubi bie ? 
anives? ragsellov ame, ons, ‘no: paeed 285) eborgam’ aber: 
fLehwitsw sty vo beqoddeeb 26w doksty bodtiom bay. \ 
naeh toma, ox d / vd) penite: eid: iva ‘Hens, epitehdate ta 

tid no bit} bao +LAbnis | /eruieaaaed uted \onigeedt ho Ys 4 


Vi rh i 


onites+ sotool fen (TER SARNEA Sig ls aupwornd: worms 
[ines # ; ; ; 1A le h - 

yie erelfuol= qebau shop 27) pilgeea ‘ede sevesot Choma 
reg Le! t I im 

ny 4 { : 4 Pa 
t6loe to Iwo veded olmanyh ats to evijsedenaey 
aie Bi : i 

soomyorsvers Deeae 


axtg' ete bamgsdiada-c fom a ty ybu as, pkg: 
ed? pn ktue ce Satiaed TiS av ispA oes ae, 


406 300 Yen bas UoOnteroe7 2 aliad FOR ety Baca 


uu 
> 
4 
- 
2, 


f 
ay 


| prises Sat Megea pnt yeen Seer, dtm o#" wer 
VORIMGSS, SIS mMATe Rey! oa Ao, epi {prone pe ani adid 
is ebonpeul Joavaluntad Por Aostamareay itdick nbd ae 
Ge bus tii al HOR oakiuedirsng, yaice baat pais | 6 
« AQaqedy as baineeey ars: eysey pnigsad teloa testucuell 
| MuCoreIatant ME ast PY Tak 
+ ie Lue ane batnoneeg b  ekeed YAdsnombns yl hob ie oni 
Laid Ri Babul ons ote a apizeisstogaade nid spanogey eo 


ne ae Sons redie’ att a5 beep 


rc Warren mates, i Saaoamer sis, sfiats tie ot 
ae “ep ‘aeesodl .a92R0L ail i spite: sz mre 


- mt r 
Are | | | 
p 7 pres ; ‘euulhiie 0 


o _ ; =e- f 1 : *, iii : 


‘a ve) Weiianes Ss 3 49 2 
Se ey en res 5 


ah 


storage bin and ducts go into the module interior and 
usually contribute to the heating requirement during the 
heating season. These are called the useful energy losses. 
However, the same losses can also cause module overheating 
during the summer and add to the cooling load. Therefore the 
present work also estimates the energy losses transferred 
from the bin walls and ducts to the space. These losses are 
presented in Chapter 4. 

In Chapter 5, the results of solar heating contribution 
assessed by different methods, collector and storage bin 
performance, dynamic behavior of storage bin and energy 
losses are discussed. This is followed by the conclusions 


drawn from this study. 
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2. The Solar Module 


2.) Introduction 

The thermal performance of the solar air heating system 
is estimated or predicted from the gathered data for module 
5 and 6 on a continuous basis. The construction details for 
modules 5 and 6, the components of the solar air heating 
System, the system operation modes, data measurement and the 


data acquisition system are presented in this chapter. 


2.2 Construction Details For Module 5 And Module 6 

Modules 5 and 6 of the test facility are idential in 
contruction except for the ceiling insulation level of 
module 6 which is about 63 percent higher than that of 
module 5. The main floor plan for each module and wall 
section are shown in Figure 2.1 and Figure 2.2 , and the 
design details are given in Table 2.1 . The construction of 
the modules represented pre 1970 construction standards for 
housing in Alberta. By uSing identical modules, module 5 
could be used as a control module to study the performance 


of a solar heating system installed in module 6. 


2.3 The Solar Heating System 

A schematic diagram of the solar ain heatingusystem 
installed in module 6 for space heating during the summer of 
1960 ise Shownein oagunes2.3. The principalycomponents of 


the solar.air-heating system are the collector, the heat 
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Figure 2.1 Main Floor Plan 
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Figure 2.2 Wall Section 
AW..B Cy ue, ae, F illustrated in Table 2.1 
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Table 2.1 Design Details Of The Module 5 And 6 


RSI - (sq-m-K/Watt) 


Floor Area(exterior) 
(interior) 


Som. Xo som 
6, 5m. 7... im 


Main Floor Wall Height : 2.4m 


Basement Height ; 


Roof (A)s 
Ceiling WB i 
Wall : 


@oénstruction (Cc): 


Wall Area ./ FOOT 


Main Floor (D): 


Basement Wall(E): 


Floor(F): 


Windows North 
South 
East : 
West 


2.4m, 1.8m below grade 


CMHC approved trusses with 0.6m-15cm stub 
210# asphalt shingles 
1.0cm plywood Ext,GD sheathing 


Standard truss with 0.6m bobtail rafters 

0.6m on center 

fiberglass insulation (RSI 5.46) module 6 
(RSet 2214 )omodule 5 

100 micrometre (4 mil) poly vapour barrier 

1.3cm drywall painted 


1.0cm Prestain Rough Tex Plywood 

8.9cm fiberglass insulation (RSI 1.76) 
Siem xae10-2em studs, 4046cm on center 

100 micrometre (4 mil) ploy vapour barrier 
1.3cm drywall painted 


Area : 1.39 


1,.6cm fireplywood 
5. lcm x] 2023cm studs, 40.6em on center 


1.3cm preservative treated plywood 

to 0.6m below grade 

5. lem eicacsinsudation (RS1vin7 6p 

to 0.6m below grade 

20.3¢m conereteswakt 

10.2cm concrete slab on 150 micrometre 
(6 mil) poly vapour barrier 


90cm x 193cm double glazed sealed 
None 

102cm x 193cm double slider 

102cm x 193cm double slider 


Window Areary “Floor -Arcast--1-.9-% 


Door 


Blectric, Duct 
Heater Capacity 


91.4cm-x 203.2cm insulated metal 


7.5 kW 
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Storage facility, the air handling system and control 
system. 

c. ihe hime ySOuEIP tecsnd solar collector tilted at 70 
degrees to the horizontal was installed adjacent to the 
south wall. The collector was assembled using 6 Solaron 
Series 2000 flat plate air collectors by connecting two in 
series and mounting three of these units in parallel. This 
assembly provided a total air heating path length of about 
3.5 m. The collector is double glazed and has an enamelled 
steel absorber. 

A 2 m’° washed gravel heat storage bin is located in the 
basement. The bin 1s of wood frame construction with 89 mm 
of fiberglass insulation and the interior surface of the bin 
has 50.8 mm of styrofoam insulation applied to it. The rock 
bed rests on the basement floor with 50.8 mm of styrofoam 
insulation between the lower plenum and the basement floor. 
Typical parameters of collector and storage bin are shown in 
Table 2.2 

The other components are an air handling unit, a 
control system with sensors and relays, motorized dampers 
for properly directing the flow, and a twopstege electric 
furnace back up heater. An internal by-pass with appropriate 
dampers in the handling unit permits a single blower to move 
air in all of the necessary system operation modes. The 
by-pass makes it possible for the blower to supply solar 
heated air to the storage bin and also to withdraw hot air 


(in the opposite direction) from the same end of the storage 
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Table 2.2 Solar Air-Heatering System Parameters 


Collector Orientation 


Collector Tilted 
From Horizontal 


eollector Gross Area‘, 
Collector Net Area 

Storage Volume And Materials 
Specific Heat Of Rock 


Storage Capacity Per 
Area Of Collector 


Air Flow Rate 


Electric Furance Capacity 


South Pacing 


70 degrees 
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bed. 


2.4 System Operation Modes 
The control system is designed to operate in five 
basically distinct modes of operation. A schematic of the 


air-flow loop for each mode is shown in Figure 2.4 


(1) Heat From Collectors 

If space heating iS required and the air temperature at 
collector outlet is 9 °C above the room temperature, then 
air is circulated through the collector, by-passes the 
storage unit and is sent directly to the space. If the 
collector outlet air temperature is less than 18 °C above 
room temperature, the first stage electric heater is used to 


keep the room temperature above 18 °C. 


(2) Heat Collection(Storage) 

If no space heating is required and the air temperature 
at collector outlet is 9 °C above the storage temperature, 
aise is circulated from.collectors to the istorage unit and 


returned ito the collector. 


(3) Heat From Storage 

If space heating is required and no energy is available 
from the collector, airiis circulated in the opposite 
direction through -shemstorage: unit and sent directly to ee 


room. If the air temperature off the storage unit is less 
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than 18 °C above the room temperature, the first stage 


electric heater is used. 


(4) Stage(Electrical) Heating 

If space heating is required, but no energy is 
available from either the collectors or the storage unit, 
then both the first and second stage heaters are in 


operation with air circulating through the storage unit. 


(5) System Purge 

If the collector temperature is above the maximum 
allowable limit, outside air can be circulated through the 
collector and exhausted to the ambient in order to prevent 


the collector overheating. 


2.5 Meteorological And Measured Data 


2.5.1 Data Measurement 


15 


Figure 2.5 shows the locations of temperature and solar 


radiation measurements for module 6. The data measurements 


can be summarized as: 


1. Electric usage is measured on an hourly basis by adding 


an Optical counter to the disk of a Standard power 
MebRer. 

2. Temperatures are measured uSing copper-constantan 
thermocouples with cold junction compensation provided 


by the data logger. 
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Til - Room Temperature 

T2 - Collector Inlet Air Temperature 

T3 "=eCollector, Outlet Air ‘Temperacure 

T4 - Air Temperature At The Top Of The Bin 

T5 - Air Temperature At The Bottom Of The Bin 
Rr e- Total Radiation On. Collector Surface 


Figure 2.5 Measurement Locations In Module 6 
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3. Horizontal surface radiation, diffuse radiation, and the 
solar radiation on the tilted collector surface are 


measured by Eppley pyranometers. 


2.5.2 Data Acquisition System 

Control of the data logging operation was accomplished 
using an Analog Device MACSYM II computer controlled data 
acquiSition system. The modes of operation for the collector 
System, the solar radiation, the electric usage, the air 
temperature for the collectors, storage bin and space were 
monitored by the data logger every two minutes. At the end 
of each hour, the logger provides totals and averages, and 
after a 24 hour period has passed, the data was transferred 
to magnetic tape in the University computer system. Then the 
memory was cleared and the sequence begins again. 

In addition, when solar energy collection occurs, a 
printed record of the 2 minute interval measurement was made 
to facilitate detailed analysis of the behavior of system 


components. 
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3. Methods For Assessing Solar Contribution 


Setalntroduction 
The solar contribution is a main parameter for the 
design and installation of a thermal system. Direct methods 
and indirect methods can be used to assess the solar 
contribution. The direct methods compute the solar 
contribution based on the measured solar records. The two 
methods are: 
1. The solar supply method: the solar contribution based on 
the energy collected by the collector. 
2. The solar required method: the solar contribution based 
on the meaSured energy supplied by the system 


components. 


For the indirect methods, 3 approaches have been taken: 

1. The building heat loss method: the solar contribution 
based on the calculated energy loss and the non-solar 
internal energy gains. 

2. The power records method: the solar contribution based 
on the power consumption of the test module and a 
control «modu leis 

3. The F-chart method: a computer program developed at the 
Uniivemsisey! of Wisconsin for theoretical predictions of 
the long-term solar!’ contribution. 


All these methods will be discussed in this chapter. 
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3.2 Direct Methods 

These methods compute the solar contribution from 
detailed measurements of the solar energy gained from the 
solar system. The data required for assessing the solar 
contribution includes the temperature difference across 
collector and storage bin during each of the operating 
modes, and electical power usage for the module. The two 
minute interval records of the electrical power usage and 
temperature difference data were summed at the end of each 
hour, and hourly values were used in the analysis. 

The hourly average temperature difference across 
collectors or storage bin can be computed by dividing the 
cummulative hourly temperature difference by the number of 


temperature readings recorded in each hour. The number of 


temperature readings recorded in each hour can be determined 


by dividing the number of minutes in each corresponding mode 


by two. The hourly sum of the temperature difference can be 


computed by mutiplying the hourly average temperature by the 


number of minutes recorded in each corresponding mode. The 


unit of time in minutes is then converted to seconds to give 


the hourly temperature differences in degree celcius second. 
The hotrhywsumeof«temperaturéerdniferencetacnoss the 
collector in each operation mode with the unit in degree 


celcius second iS given by: 
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COL 
DTCcol = { L are) = 2 x60 C32 4") 
1 
HFC 
DTChic Maa x ATC pextotxeco (342) 


2 


And the hourly sum of temperature difference across storage 


bin for each operation mode is given by: 


EH 
DrSeny t= (Gee ars.) x 2) x 60 Cees 
ay 
HFS 
BTShESH= US TIATS Pour 2yx<¢60 (B24) 
qt 
COL 
DTScolyesk Vearssay x 2 x 60 Cor S)) 


al 


Where 

EH = no. of minutes in the electric heating mode 

HFS = no. of minutes in the heat-from-storage mode 

HFC = no. of minutes in the heat-from-collector mode 

COL = no. of minutes in the heat-collection(storage) mode 

aTC = instantaneous air temperature difference measured 
across collector, Tfe-Tfi, [°c] 

dTS = instantaneous air temperature difference measured 


Bcmoecestorage mane arstersh ed °Ca 
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3.2.1 Solar Contribution Using Solar Supply Method 

The energy collected by the solar collector is the 
maximum solar energy availiable to Module 6. The collector 
Operates into 2 modes: the heat-from-collector mode in which 
the collected energy is sent directly to space heating, and 
the heat-collection(storage) mode with the collected energy 
sent directly to the storage unit. The maximum energy 
Supplied is the sum of the energy delivered from the solar 
collector in each operation mode and the electrical energy 
Supplied. The hourly collected energy can be estimated by 
multiplying the hourly sum of temperature difference across 
collectors by the product of the mass flow rate and fluid 
heat capacity. The maximum hourly energy supplied to space 
heating can be estimated by: 


Qhsun Ochfc +.Qccol (3.6) 


Qhsun + Qeh wee yes 


iH} 


Qht 


Where 
Qhsun = hourly solar energy gained from the collector, [J] 
OQchfc = hourly energy gained from collector 
for space heating(heat-from-collector mode), 
menUDTehic) sui 
Qccol = hourly energy gained from collector for storage 


(heat-collection mode), mCp(DTCcol), [J] 

Qeh = hourly energy supplied from electric (stage) 
heating mode, [J] 

Qht = hourly total energy supplied for space heating, [J] 
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The cumulative monthly solar energy supplied for space 


heating (Qmsun) can be determined from: 


K K 
Qmsun = Z Qchfe + Z Qccol (328) 
Ne k 


The cumulative monthly total electric energy usage 


(Qmeh) is: 


K 
Qmeh = £ Qeh (3:9) 
k 


Therefore the maximum solar contribution supplied to 


Space heating can be estimated as: 


Omsun 
(3.10) 


nS 
Omsun + Qmeh 


3.2.2 Solar Contribution Using Solar Required Method 

Using this method, the solar contribution can be 
assessed by computing the total energy required for space 
heating by an energy balance for the module. Total energy 
required for space heating can be estimated by summing the 
collected energy from the collector that is immediately used 
for space heating, the energy supplied from storage bin, the 
energy supplied from storage bin during electrical heating 


and the electrical energy usage. Thus, the hourly total 
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energy reguired for space heating can be written as: 


Qht = Qchfc + Qshfs + Qseh + Qeh (B11) 
Where 
Qht = hourly total energy required for space heating, [J] 
Qshfs= hourly energy supplied from storage bin for space 
heating (heat-from-storage mode), mCp(DTShfs), [J] 
Qseh = hourly energy supplied from storage bin with 


auxiliary heater operating (heat from storage bin 
during electric heating mode), mCp(DTSeh), [J] 


Ochfc and Qeh are the same as declared in Equations 3.6 
and soi. 


The cumulative monthly total energy required for space 


heating can be determined by: 


K 
Ome =) 2e0Ohe C3512) 
k 


By combining Equations 3.9 and 3.12 , the monthly solar 
contribution required for space heating can be assessed as 


follows: 
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3.3 Indirect Methods 


3.3.1 Solar Contribution Using Building Heat Loss Method 

In the process of transferring solar energy from the 
collectors, significant energy losses may occur for active 
solar systems. As an alternative to using detailed 
measurements to determine the energy losses of an active air 
solar system, a method studied by Taylor[16] called the 
building heat loss method can be modified to evaluate the 
solar space heating contribution of the solar system. This 
method, which defines the building(module) as the analysis 
control volume, can be used to assess the solar contribution 
with a minimum of instrumentation and data acquisition. 

The overall solar contribution for space heating can be 


computed as follows: 


Overall Building Non-solar 
solar = heat + internal C3212) 
contribution loss energy gains 


The overall! solar ‘contribution 15 thegoum or sche ner 
passive energy gains (Qsp) and the solar energy supplied to 
the building interior by the active system (Qsa). Direct 
passive windows gains of small magnitude can be estimated by 
using measured radiation data. The estimated passive gain is 
substracted from the computed overall solar contribution to 


obtain the active system contribution. 
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Non-solar internal energy gain is the electrical power 
supplied to the module which includes: 

1. The power for electrical heater 
z-unTherpoWernforylight ingpeand 
3. The power for the heating system fan. 

Since these three items contribute to space heating, 
the non-solar internal energy gain was estimated to be the 
total electrical energy recorded by the power meter. 

The building heat loss can be estimated from the 
building heat loss coefficient (UA) times the difference 
between the module interior temperature and the ambient air 
temperature. Therefore, the hourly energy loss for the 


module can be written as: 


Obhia=s (UA) &TEeStTA) (3615 ) 


Equation 3.14 can be expressed as: 


Osaittidspw= lObhbhm= deh (3716) 


There are no south windows in module 6 and shading of 
the other windows by adjacent buildings reduces the direct 
passive solar gain. Also, aluminum foil was added to the 
windows to further reduce this solar passive gain. Indirect 
passive solar gains due to solar radiation absorbed by the 


walls and roof are difficult to estimate as the thermal 
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performance of the components is altered. At any rate, 
passive gains will result in a reduced conduction energy 
transfer to the outdoors, or, for large passive gains, a 
transfer of energy to the module. For module 6, the attic 
Space is well ventilated and the solar collector shades most 
of the south wall. In addition, a different indirect passive 
solar gain occurs for each month. Calculation of the effect 
of the indirect passive solar gain was not possible and its 
effect on the heating load was assumed to be zero. Then 


Equation 3.16 can be rewritten as: 


Qsa = Qbhl - Qeh (32172) 


To minimize the effect of the indirect pasSive gains, 
the building heat loss coefficient (UA factor) in Equation 
3.15 was calculated using data in December and January for 
wWhichealitt le or no bright; sunshine hours were recorded. In 
these months, solar energy iS a minimum and the energy for 
heating the module will come from the heating utility . The 
UA factor 1s estimated from Equation 3.17 bysneglecting the 


solar term and the UA equation can be written as; 
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If the whole month of December or January cannot be 
used as the analysis period, days in each month with little 
or no solar energy can be used to estimate the UA factor. 


Then the UA equation can be written as: 


M 
UA = { L —————— } / M Gor 19) 
m 


Where M is number of days in the analysis month with 
little or no solar energy 
The monthly total active solar gains (Qmsa) can be 


computed by: 


K 
Onsae= 2 Obhls 2: Gen (3920) 
k 


And the fraction of total energy for space heating 


Supplied by the active air thermal system is given by: 
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K K 
Qmsa B- ODnL =" Oeh 
k k 
nbhl ERIE Oo SRE (mon 
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Brn 2 Solar Contribution Using Power Records Method 
Comparing the energy usage by two buildings, module 5 
and module 6, is another method for expressing the 
contribution to space heating by the solar energy system. 
These two modules have identical construction except for the 
RSI value of the ceiling insulation level of module 6 which 
is 63 % greater than that of module 5. Thus, a comparison of 
electrical energy requirements for each module should 
indicate the solar contribution for space heating. The power 


ratio Py of modules 5 and 6 is given by: 


K 
Cumulative monthly electric Z Qeh6 
power usSage in module 6 k 
Py = —A —_ = om2 ce 
Cumulative monthly electric K 
power usage in module 5 x Qeh5 
k 


Where 
Qeh5 = hourly electric power usage in module 5 
Qeh6 = hourly electric power usage in module 6 


no. of hours with electric power usage in a month 
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However, although modules 5 and 6 have identical 
contruction, the UA factors for these two modules may 
differ. The power ratio Py of modules 5 and 6 is affected by 
any deviation of these UA factors and may be compensated for 
by calculating the power correction factor of these two 
modules, Fua, under conditions with little or no sunshine. 
Reom the discussionn,of UA factor tn section 3.3.4) ) “the 
monthly power correction factor can be estimated by using 
Equation 3.19 for each module. 

Monthly UA factor of module 6 


Fua = — = ec) 
Monthly, UA factor of module 5 


By “combining Equatrons 13922 Card (352867. ,the monthly 


actual power ratio can be found as: 


Pya = Py / Fua 


Therefore the monthly fraction of energy for space 


heating supplied by the solar system can be found as: 


meh = 1 = Pya (ame } 


3.3.3 Solar Contribution Using F-chart ‘Method 
In the development of the computer program called 
F-chart, assumptions were made that the systems are 


well-built, flow distribution.to collectors is uniform, flow 
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rates are assumed constant, and the system configurations 
will be close to those for which the correlations were 
developed[27]. This computer program can be used to predict 
the long-term solar heating system performance. 


Two dimensionless groups, X and Y, were identified.[27] 


Ref.collectos.loss AcFrUL(Tref - Ta)Ac 
a (3225) 
Heating Load L 


Absorbed solar energy AcFr(ra)eHt.N 


Aes G2s26) 
Heating Load L 

Where 
Ac = collector area, [m?] 
Fr = collector-heat exchanger efficiency factor 
UL = collector overall loss coefficient, 

Wyn 2G 
At = total no. of seconds in the month 
Ta = monthly average ambient temperature, [°C] 
Tref = an empirically derived ref. temp., [100 °C] 
L = monthly total heating load for space heating, [J] 
Ht = monthly average daily radiation incident on 

the collector surface per unit area, [J/m?] 

= days in month 

(Ta)e = monthly average transmittance-absorptance product 


A heat exchanger is not required for an air system. The 
term Fr is replaced by the solar collector heat removal 
factor Fr, and if the duct losses are assumed very small, 


Equations 3.25 ende3.26 can be rewritten as: 
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Frigadvec xi ft .N © Ac '/ L (3.226) 


A dimensionless factor f, called the monthly fraction 
of heating loads supplied by the solar system, relates these 


two dimensionless groups X and Y is expressed as follows: 


m= EROn yo 0706529> 071 59Y" 220700 187 xe re 000Sy= Co29) 


Theyr-chart predictions of solar contribution for the 
system installed in module 6 were based on the 
meteorological data obtained from the International Airport 
and the long-term solar radiation records obtained from 


Stony Plain about 60 km west of Edmonton. 
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4, Rating The Thermal System 


4.1 Introduction 

Many designers and potential users of solar heating 
Systems developed different methods of predicting and 
measuring the performance of the thermal system. The 
collector and heat storage facility are the principal 
components in a Solar air-heating system, therefore, their 
rating is important. Rating can be broadly categorized as: 
1. Thermal performance 
ase COSt 
3. Life time, durability, maintance and ease of 

installation. 

For rating solar air collectors, Kenna[38] pointed out 
that the rating should take into account the economics and 
durability factor , together with the thermal performance. 
The thermal performance of a collector 1s compared by the 
concept of collector thermal efficiency and has been studied 
by a number of investigators, Hottel and Woertz[39], 
Bliss(40)> Whiblirer 40) and Du€éiie and Beekman! 26) (27) (tial) 
and "Streed[33], Lof; Shaw and Oank[17], Siebers and 
Viskantai42], and Wunde[8) ‘calculated the tpertormance ‘of 
solar collector based on meteorological data. A number of 
reseachers, Klein, Cooper, Freeman, Beckman and Duffie[12], 
predicted the performance by various theoretical models. 

For rating storage beds, the measured and predicted 


performance of rock bed storage has been studied by Persons, 
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Duffie and Mitchell[18], Hughes, Klein and Close[19]. A 
mathematical model was developed for the rock bed and used 
to study the performance of solar heating systems and 
cooling systems. 

The overall performance of a solar heating system is 
the result of the thermal performance of the collector and 
storage bed in conjunction with the control system used. The 
evaluation of system and component performance will be 


discussed in this chapter. 


4.2 Rating The Solar Air Collector By Thermal Efficiency 
The steady state thermal performance of the solar 

collectors is determined by measuring the energy collected 
for a combination of values of incident radiation, ambient 
temperature, and inlet fluid temperature. The performance. 
can be described by an energy balance that indicates the 
distribution of incident solar energy into useful energy 
gain and various losses. The energy balance on the collector 


can be written as: 


The rate of The rate of The rate of The rate 
useful energy energy loss from of energy 
- energy = absorbed in = ‘the collector - storage 
extracted the collector by conduction, in 
from the by the convection and the 
collector absorber radiation collector 


The rate of useful energy extracted from the collector 


is equal to the rate of energy carried away by the transfer 
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fluid. If the collector is at steady state, the rate of 


change of energy storage in the collector materials is 


assumed to be zero, and the energy balance equation becomes: 


Qu m 
Hitcoveestuil( Gp) =bhtabosisetpt(Tre - TE) (4.1) 
Ac 


To preclude the necessity for determining the average 
it 1S convenient to introduce 


temperature of the absorber, 


the parameter Fr, where 


actual useful energy collected by 
a flat=pilate collector 


Fr = 
useful energy collected if the entire flat-plate 
collector surface were at the inlet fluid temperature 


into Equation(4.1) results in 


InteoducmnghFry 
Qu m 

Sehr LHe ra)e —LU CRE seta) =3-—Cpo (tie = <Thi) (4.2) 
Ac Ac 


This equation is called Hottel-Whillier-Bliss equation 


(Simth[20]). It expresses heat delivery from a collector in 


terms of two operating variables: 
diffuse and reflected) radiation 


The total (beam, 
and 


1s 

incidenewom the collector cover plate, Ht, 
2. The temperature between the entering heat removal fluid, 
and the ambient temperature, Ta. 
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4.3 Governing Thermal Efficiency Equation 


The solar collector efficiency is defined as: 


Actual useful energy collected Qu/Ac 
fi = 


Solar energy incident upon or Ht 
intercepted by the collector 


According to Equation(4.2), the collector instantaneous 


efficiency of the flat-plate collector is given by: 


GRe hte a) mep (Tier wia) 
fi = Fr[(ra)e - UL ——————] = ————_____ (4.3) 
Ht Ac.Ht 


Equation(4.3) indicates that if Fr(ra)e and FrU! are 
eonstant, the plot of 61 versus (Tf£i-Ta)/7Ht. is a Straight 
line with the slope equal to FrU! and the y-intercept is 


equal to Fr(ra)e. 


4.4 Computation Of Collector Thermal Efficiency 

The collector efficiency discussed in Section 4.2 and 
4.3 is the steady state instantaneous efficiency. However 
instantaneous readings are impractical because actual 
thermal systems operate under variable conditions. The 
thermal efficiency is very sensitive to a sudden change of 
solar radiation due to cloud, ambient temperature and wind 
velocity. This suggests that the time averaged collector 
efficiency over selected periods of time be used in 


calculations to compensate for the variation in operating 
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conditions. The collector efficiency is often presented as: 


(A) Collector efficiency on an hourly basis 
Hourly collector efficiency is selected because it is 
the minimum time period for meteorological records. The 
efficiency of the collector during each hour of operation 
can be found from 
Ou mep(T£e = Tri) 


gbirs Sl Seer sah- Soest (4.4) 
Ht .Ac Ht.Ac 


Where Ht is the radiation on the collector titled 
surface , and Qu is the rate of useful energy gain for the 


ROUr, 


(B) Collector efficiency on a daily basis 
The daily efficiency, studied by Duffie and 
Beckman[26], is not the average hourly efficiency but must 


be calculated from 


I iL 
Ou mep 2 (rre = Thi) 
t i 
£d5= = ao) 
I I 
AChavHs ACeZ ant 
i z 


Where 2 Ou ¢s the total useful energy ofor the 
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I 
day and = Ht is the total incident solar radiation on 
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the tilted surface 


(C) Collector efficiency on a monthy basis 

Both the hourly and daily efficiency may still be too 
short a time period to not be affected by the instantaneous 
meterological conditions. Therefore the monthly collector 
efficiency has been suggested as the best assessment of the 
thermal performance of the collector for system design. The 


monthy efficiency can be calculated by 


K K 
x Ou ncCpee (The + TLL) 
k k 
En = (4.6) 
K K 
AG 1DGHt Acez Ht 
k k 
K 
Where £ Qu is the total useful energy for the 
k 
K 
month and £ Ht is the total incident solar radiation on 
k 


the tilted surface. 


4.5 The Performance Of The Storage Bin 

Alihought neat from solar air collectors cam be (stored 
in tanks) Of water or other liquids, in phase-change 
compounds, in numerous small containers of liquids, and in 
bins of loose solids, only a pebble bed is technically and. 
economically advantageous[15]. The storage bin is one of the 


principal components in the solar air-heating system and the 
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performance of the storage bin will directly affect the 

performance of the whole thermal system. Therefore, rating 

the performance of the storage bin is important. 

The storage unite for this thermal system is a rock bed 
with 2m* of washed gravel enclosed in an insulated wood 
frame box. The energy storage capacity of the bin per area 
Giecollector is 545.6 kj/m’—="C } and thesbin can store all 
of the solar heat collectible on a cloudless day. The 
storage bin operates in three modes: 

1. Collected energy goes directly to the storage bin. This 
is called heat-collection mode which is symbolized by 
ECO. 

2. Stored energy from the bin directly supplies the space 
heating, this is called heat-from-storage mode which is 
Symbolized by ‘hfs'. 

3. Stored energy from the bin and the electric heater 
supplies the space heating. This mode is symbolized by 
Nene, 

Figure 4.1 shows the schematic diagram for the charging 
and discharging direction of the storage bin. During 
charging, the air flows in the downward direction with the 
aire temperatunewat the top-of bain, (Tst)yhigher= than the aim 
temperature at the bottom of bin (Tsb). During discharging, 
thewairetlowss2a the Upward direction, 

From the theoretical conservation of energy point of 


view, the energy balance on the whole rock bed can be 


written as: 
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Heat From Heat 
Storage Mode Collection Mode 
(Discharging) (Charging) 
a ae Tst 
Storage Bed 
Dae cer Tsb 


Heat From Storage 
During Electric 
Heating Mode 
(Discharging) 


Tst - Air Temperature At The Top Of The Storage Bin 


Tsb - Air Temperature At The Bottom Of The Storage Bin 


Figure 4.1 Direction Of Charging And Discharging Of The 
Storage Bin Corresponding To Each Operation 
Mode 
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The rate of The rate of energy The rate of 
energy leaving the bed energy storage 
entering the - (include airleakage = in the rock 


rock bed and conduction bed 
through the walls) 


In equation form: 


Qscol - (Oshfs + QOseh + QOsloss) Osto Care) 


Where 


Oscol = rate of solar energy entering the storage bin 
Eeomacollector “during heat-collection mode col”|, 

: mep(DTScol)., [Wl 

Oshfs = rate of storage energy from the bin supplied 
for space heating during heat-from-Sstorage 

; mode nis’, mep(DTshtis), [WwW] 

Qseh = rate of storage energy from the bin 
for space heating during electrical heating 
mode'eh', mCp(DTSeh), [W] 


QOsto = rate of energy accumulated in the storage 
t bin, [Ww] 
Qsloss = rate of storage energy losses from storage bin 


through conduction or air leakage, [W] 


The monthly amount of energy loss and storage in the 


storage bin is: 


K K K K K 
BuOscoly-TaOshis SBELOSeEh SSLSOslossP+EasOsto (4.8) 
k k k k k 


According to the ASHRAE proposed method for the energy 


balance of the storage bin, the net energy accumulated or 


i 


16: 8Tea ene “gh a ap 54 a 
sosiors ‘Tessas bat a. poly . 
A507) SASHES Paeetape ts ie ne a ol 
fied - AMUEAOS BAB. 1), 0 

: a ada: Revoxds | 


te 
I j 
~ ‘Ry 7702 
Fy 
) / 
of20) = (meolad > Need + 2 2ae8) 7 


id somota 4d2 sntastns donans\ setoa We edad’ = Os 
‘tes ebont, norgse lle ar sa haan paint alah == sain a 
Ul ‘ P40! 
ba hice Grd. sng i) 
sye2ote-iig- ss 


hae si ites 
; i } ee en : : 
pidedd mMets Voss abetate tO ete1 = 
oni deen a ia 2 thai ett w! ah See 70% 1. oa 
{at i |e" Sbomt = 
soaTogR) etd iva ne ai Ae *e2a1 * O32 


Vad GT 7 
ss sei ede a 


. ete 


aie i a alm -. 

sasha epzi6ie hs = a ons shui oer ott Mo 

- | jai. wid opsit je 

“¥ | : ies 
(648) o22Q "i * E eid site> 2 - Isom 
in 4 ov | i 


a 


4 | 


stored inside the rock bin for a monthly period is assumed 


to be zero. Therefore, Equation 4.8 can be rewritten as: 


K K K K 
Z Qsloss = Z Qscol - Z Qshfs - Z Qseh (4.9) 
k k k k 
Then 
Monthly energy losses from the bin 
LS = 


Monthly solar energy stored in the bin 


K 
LVOSloss 


SS (4.10) 


2 OSCOL 


Equation 4.10 is the fraction of monthly energy lost 


from the bin through conduction or air leakage. 


4.6 Energy Transfer From Ducts And Storage Bin 

Energy losses or energy transfers from the ducts and 
storage bin should be determined because they will directly 
affect the performance of the thermal system. When the 
building requires heating, losses from the storage bin and 
ducts that go into the heated space are useful energy 
losses. But losses from the duct section located outside the 


heated space are detrimental energy losses. In order to 
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estimate the energy transfers from the thermal system to the 
module, an energy balance for the module and the thermal 
system is needed. The monthly solar energy gained from 
collectors and the monthly total energy required for space 
heating have been discussed in Chapter 3. The monthly actual 
solar energy supplied from the thermal system to space 
heating, Qmas, can be determined by substracting the monthly 
electrical energy usage from the monthly total energy 


required for space heating. 

Qma's = Qmtr- Qmeh Catt) 

Then the monthly energy transfers or losses from the 
thermal system, Qsysl, can be determined by subtracting the 
monthly actual solar energy (Qmas) from the monthly total 
solar energy (Qmsun) gained from the collector. 

Qsysl = Qmsun - Qmas (4.12) 
and 


Lsys = Qsysl / Qmsun (Apaes:) 


Equat jone4y 1Seisotheifractioneofivenergy i transfterseor 


losses from the thermal system in each month. 
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5. Results And Discussion 


wb antroduction 

The characteristics of the solar air heating system 
were investigated using the methods developed in Chapter 3 
and 4. The solar contribution to the heating load was 
assessed byidirect and andirect methods, the thermal 
efficiency of the collector calculated for different time 
periods, and the dynamic behavior and the thermal 
performance of the storage bin examined. System ductwork and 


storage bin energy transfers will also discussed. 


5.2 Direct Solar Method 

The direct solar method is based on measured energy 
quantities to assess the solar contribution. Therefore, the 
monthly total number of hours in each system operation mode 
is presented in Table 5.1 . Total hours of collector 
operation is the sum of hours in heat-from-collector mode 
and hours in heat-storage mode. The results indicate that 
the bed temperatures are too low for operation in 
heat-from-storage mode during December 1981 and January 
VOS25. 

The direct solar methods can be divided into 2 parts, 
the solar supply method which is based on all the collected 
energy supplied by the collector and the solar required 
method which is based on the energy required by the module. 


Therefore , energy losses are considered in the latter. The 
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Table 5.1 
Month Yr 
Sep 8 1 
Oct 8 J 
Nov 8 1 
Dec 8 1 
Jan 82 
Feb 82 
Mar 82 


Season 


Number Of Hours In Each System Operation Mode 


Heat 


From 
Collector 


Mode 


HEE 


(Hr ) 
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Storage 
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fraction of energy supplied by the solar system was 
calculated on a monthly basis and for the heating season. 
The results are presented in Tables 5.2 and 5.3, and 
indicate that 11.8 to’ 15.7 % of ‘total energy was ‘supplied by 
the thermal solar system for the heating season from 
September 1981 to March 1982 . The energy supplied to the 
module is supposedly equal to the energy required by the 
module. The unequal result for these two methods may be due 
to the energy losses and energy transfer from ducts and the 
storage bin. 

The solar contribution in the month of September 1981, 
62.8 %, was obtained from the supply method while 29.3 % was 
obtained from the required method. There is a large 
difference between these two monthly solar contributions. It 
appears to be due to the fact that September 1981 was a very 
warm month. The collected energy directly sent to space 
heating is only about 14.0 % of the energy sent to the 
storage unit. Also, the energy required from the storage 
Unit for Space heating was only about one half of the total 
energy stored in the bin. The temperature inside the bin was 
very (high andyresultéd in an increase of energy transter 
through the walls of storage bin to the basement. Thus, the 
large discrepancy between the solar supply method and the 
solar required method is mainly caused by energy loss from 


the solar system ductwork and loss from the storage bin. 
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Table 5.2 Monthly Solar Contribution Assessed By 
Solar Supply Method 


Month Yr Monthly Monthly Monthly Percent 
Electric Incident Heating Solar 
Energy Solar Load 
Qmeh Omsun Omt nS 
(MJ) (MJ) (MJ) (%) 
Sep 8 1 T6305 aes 0 3104.3 625.0 
Oct 8 1 3927, 5 904.4 AGSAi So MERAY 
Nov 8 1 BOOS 9027.3 6673.6 ars 9 
Dec 8 1 Oi. 2 S70 10404.2 8.4 
Jan 82 1326455 669.3 U3 235 40 4.8 
Feb 82 95 UG 1475.0 LO9s2 02 1376 
Mar 82 6418.7 2129.2 854739 26.0 


Season 58478.9 evs i 
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Table 5.3 


Month Yr 
Sep 81 
Oct 8 1 
Nov 8 1 
Dec 8 1 
Jan 82 
Feb 82 
Mar 82 


Season 


Monthly Solar Contribution Assessed By 


Solar Required Method 


Monthly 
Electric 
Energy 


Omeh 
(MJ) 
16.21% 
3927. 
5691 
9533. 
13264. 
9508. 


6418. 
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5.3 Indirect Solar Method 


5.3.1 Building Heat Loss Method 

The building heat loss method is an alternative method 
of assessing the solar contribution by using a control 
volume analysis of the module. The advantage of this method 
is that the solar contribution can be assessed using a 
minimum of instrumentation and meteorological data. The 
parameters used are: the building heat loss coefficient of 
the module, the room and ambient temperature, and the 
electrical energy usage. The building heat loss coefficient 
was calculated to be about 130 W/°C in the month of January 
1982, which is the month with little solar energy supplied 
from the solar thermal system. The solar contribution was 
assessed by substracting the energy losses from the module 
from the electrical energy usage and the results are 
presented in Table 5.4 . During the heating season, 13.8 % 
of total energy was supplied by the solar system and the 
monthly solar contribution ranged from 040 48.2%. 4Thevzeno 
contribution occurred in January 1982 and appears to be due 
to damper leakage since the measured solar contribution for 
January was about 4.6 %. The night time hourly temperature 
records show that the air temperature at the collector inlet 
is always greater than the air temperature at the collector 
outlet. This indicates that leakage is occurring but system 


configuration did not allow measurement of this leakage. 
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Table 5.4 


Month Yr 
Sep 81 
Oct 8 4 
Nov 81 
Dec 8 4 
Jan 82 
Feb 82 
Mar 82 
Season 


Monthly Solar Contribution Assessed By 
Building Heat Loss Method 


Monthly 
Electric 


Monthly 
Heating 
Load 
Qbhl 
(MJ) 
3148.6 
5519. 
6804.0 
10332.4 
£32765. 0 
10694.5 


6239.4 


57964.0 


Percent 
Solar 


nbhl 
(%) 


48.2 
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A low leakage damper on the outlet side of the solar 
collector isolates the collector when solar energy is 
unavailable. Fan operation during the heating cycle places 
the collectors under pressure and any leakage past the 
damper introduces cold air into the module. This 
characteristic is inherent in the system design and its 
effect must be kept to a minimum by good damper seals and 


careful adjustment of the control linkage. 


5.3.2 Power Records Method 

The power records method assesses the solar 
contribution by comparing the monthly power usage for module 
5 and 6, and then calculating the energy supplied by the 
solar system. The heat losses for module 5 and 6 are not 
identical, but a correction factor 0.97 (Appendix A-2), was 
added to the power equation to adjust for their deviation. 
The monthly solar contribution is presented in Table 5.5 
The table indicates that 11.3 % of total energy was 
contributed by the solar system during the heating season, 
and the monthly energy saving ranged from 0 to 43.1 % . The 


January solar contribution was zero. 


5.3.3 F-chart Method 

F-chart is a versatile, self-explanatory interactive 
computer program developed at the University of Wisconsin 
for the analysis and design of active solar heating systems. 


Using the parameters of the thermal systems, the 
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Table 5.5 


Month Yr 
Sep 8 1 
Ot 8 1 
Nov 8 1 
Dec 8 1 
Jan 82 
Feb 82 
Mar 82 


Season 


Monthly Solar Contribution Assessed By 


Power Records Method 


Monthly 
Electric 
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meteorological data from the International Airport and solar 
radiation data from Stony Plain near Edmonton in the 
computer program, the long-term average system performance 
was calculated and is presented in Table 5.6. The results 
indicate that 17.8 % of total energy was contributed by the 
solar system in the heating season, and the monthly solar 


contribution ranged from 10 to 44 % 


5.4 Summary Of Solar Contribution Methods 

The monthly solar contribution assessed by different 
methods for the heating season from Sepember 1981 to March 
1982 is presented in Table 5.7 . The results indicate that 
the contribution varies from 11.3 % assessed by the power 
records method to the 17.8 % contribution predicted by the 
F-chart using the meteorological data from the International 
Airport and Stony Plain. In Table 5.7 , the solar 
contribution 15.7 -% calculated from the solar supply method 
is higher than the 11.8 % obtained from the solar required 
method. That is an expected result because the deviation of 
about 3.9 % may be due to the energy losses of the thermal 
system. As mentioned in section 5.2 , the contribution of 
September 1981 calculated from the supply method is about 
62.8 % which ismuch higher than 29.3 % obtained from the 
required method. The main reason for the difference appears 
to be the energy tranfers from duct work and the storage 
bin. The data undiestes that only 14 % of monthly Solar 


energy collected was used immediately for space heating, 


taloe bie J10qziA thaeidonisisnit any mn7? agso 
eid mi noscombe 1aen nisl YNosrk mew eset 
arnamiotreq neteye spsis¥e mIst~onD fl ‘enz 2 IRE 
adfyes: oAT . 8.8 eldeP ab bettsessq ef bas & 


siz ud betudtiiaeo @aw Yves tates 30 x & cf Sede! ot 


xefoe yindaom sit Bne faoased enissed a2 ns me 2eye. 
: : ‘ . nO 
° ws 52°01 mori Saphss nocsediae 


~ 


shodyeM aoltgdisigaed pee LeP, x6 Se leat 


inerethib vd Bearezes ani sudingeos ‘Bios y teams ale 
‘(owaM of TSG! vedmegee mort, “bebse preteen ads) 262 
‘ov atedsiba? acluest) an? . (2 sidet in? hetnessig at 
~awoq edt yd hoeseeee 2 BE. 1 i) mini? aomey heii seed Sa 
siz yd betotSexq ‘oPsddi vanes £5. % sds’ os cer am 
[snoiisaseinl ot, woTt aren te spo loudetem Sci 5 i 
sgsloe of? « Vad SigaT om! tele vans2d* his 3 > 
‘bodtem giggue isto? Sms mort betelysiss @ 5,€) mossemegs 
betiupe? xeaice edo mo tl Haniesdo Pa a ais reddy 2ate 
26 aottstveb eit saueced Giwsax bavAbasdila: £2 rartT es 
{gmient enia te esaeol ypiane Bis. Ot SEO ad Lei > t 
to. noisudiseno> efi , $36 Aoliose on baat lamem 24 sia 
suods ai Sod dan viggie, ai me's? Otisluoise Se! ted 
ad? mov? hantagao 2 €528 mais wedeit dau et doidy 2B. 
steeeue sane 7etiit ede yo? nnasey Aiem.sdT .boiten rr 
ogaeigse S17 bas /7oM, Fob Mar? eteinsit yptses sda. ei 
weloe qidinon to & &) qinc dedt eezesthit, ered al 
Gas oaen soage 10% isreibenm: beau sw aaa 


a 


a, 


SS 


Table 5.6 Monthly Solar Contribution Assessed By 
F-chart Method 


Month Yr Monthly Monthly Percent 
Incident Heating Solar 

Solar Load f 

(MJ) (MJ) (%) 
Sep 8 1 4920.0 2800.0 44.0 
Oct 8 1 4840.0 4820.0 Se8 
Nov 8 1 SogUn U 7940.0 Or 
Dec 8 1 3450.0 10680.0 T0330 
Jan 24 Sonn. 0 ibe AOC, 20 
Feb 82 4910.0 9580.0 18 
Mar 82 6300.0 SE70).0 22h) 
Season 56660.0 PANS) 


The #@Parameters Used: Fruh = 2.9 W/m4-4C 
Fri¢rage = 0/25 
UASFactor =)150 20 W/ 7G 

Collector Flow Rate x Specific Heat / Area = 13.95 W/m?-°C 


Storage Capacity Per Area of Collector = 545.6 kJ/m?-°C 
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Table 5.7 Monthly Solar Contribution Assessed By 
Different Methods 


Month Yr Direct Methods Indirect Methods 
Percent Percent Percent Percent Percent 
Solar Solar Solar Solar Solar 
ns nx nbhl neh f£ 
(%) 6%.) G%) &%0) ey) 
Sep 8 1 6248 29 = 3 48.2 4321 44.0 
Oct 8 1 1942 L625 28.8 30 «2 3 ta0 
Nov 8 1 14.9 17 16.4 hee0 Waste 
Dec 8 1 8.4 2640 cea” aa2 #O:,:0 
Jan 82 4:8 4.6 O&O O20 h280 
Feb 82 1356 (2 su hdd Sear hh <0 
Mar 82 26.8 nO .22 2281 22°6 2280 
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while the rest went to the storage unit. However the energy 
losses through ducts and walls of the storage unit are very 
Significant and this will be discussed in section 5.9 . It 
shows that the more solar energy collected in the storage 
unit, the greater is the deviation between the two methods. 

For the indirect methods, the building heat loss method 
and the power records method are very similar. Both of these 
methods can be used to estimate the solar contribution 
without using collector operation and radiation data. The 
difference between these two methods is that the building 
heat loss method uses the test module as the control module 
while the power records method requires a second module 
having Similar thermal characteristics. The main parameter 
of these two methods is the overall building heat loss 
coefficient which is difficult to determine but is assumed 
to be constant for a specified period. The results obtained 
from these two methods for the heating season indicate that 
the solar contribution is only 0.5 to 2 percentage points 
different from that of solar required method. The solar 
energy collected was calculated from records of collector 
operation and should be reliable. The January contribution 
of about 5 % is reasonable as 60 hours of collection 
OCCUrred during thisimonth. It 1ssalso acceptable) that soler 
contribution calculated by the building heat loss and the 
power records method was zero during January. Damper leakage 
during a cloudy cold January would increase the power 


consumption. 
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F-chart is a theoretical prediction of the long-term 
solar Tcontributroen. [Thetsobartcontrirbut ron sofclyi78 cy incthe 
heating season obtained from this method is higher than that 
obtained from the other methods. Both of the F-chart 
predictions and the solar supply method assumed no energy 
losses in the thermal system. 

The monthly solar contribution for the heating season 
indicates that: 

1. The indirect methods are acceptable methods to assess 
the solar contribution. 

2. Long-term monthly solar assessment by the indirect 
methods appears to be closer to the actual measured 
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5.5 Performance Of The Solar Collector 

The published performance curves for a solar collector 
represent its performance under steady state conditions. But 
solar systems operate under dynamic conditions and the 
equalities defined by Equations 4.1 and 4.2 do not govern 
the thermal performance of the collector since part of the 
solar energy absorbed is used to heat the collector and its 
components. Because the actual collector performance is 
under dynamic conditions rather than under steady state 
conditions, comparison between the actual performance and 


the steady state performance 1s necessary. 


5.5.1 ASHRAE Standard Test Method 

The Solaron 2000 solar collector used on module 6 had 
been rated using a test method similar to the ASHRAE 
proposed standard test. But the test data were presented as 
a plot of collector efficiency versus (Tfe-Ta)/Ht, where Tfe 
is the air temperature at the collector outlet. Conversions 
Were) applied Loseonvertythe test result to theystandard 
parameters referenced by Duffie and Beckman[27], Fr(ra), and 
Frul, this gave: 
1, At nommaiwoncidencesrrizva), = Er(ra)< =—0s5 
on) Bheveollectormieatwsuoss. taco, Ew ee Wl se oC 

For anwarp atest flow rate per unit erea of collector of 
0.0102 m?/s-m?, the resulting efficiency curve is shown in 
Figure 5.1. With’ tiwe=theoretical performance.of the 


collector known, one may ask, what is the actual performance 
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of the Module 6 collector constructed using these units and 


operated in a northern climate? 


5.5.2 Instantaneous Efficiency 

The instantaneous efficiency of the solar air collector 
under a cloudless sky was studied by selecting two days in 
March 1982. The collector operated dynamically and the 
collector efficiency was calculated for each two minute 
interval record over the operation period for that day. The 
results are presented in two plots. Figures 5.2 and 5.4 
shows collector efficiency versus (Tfi-Ta)/Ht. Under steady 
State conditions, a straight line should result. Figures: 5.73 
and 5.5 present the solar radiation on the collector tilted 
Surface and collector efficiency versus the time of day, and 
illustrate the trends of solar radiation and collector 
efficiency during the collection period. The number on the 
horizontal axis indicates the beginning of the hour. 

The efficiency curve in Figures 5.2 and 5.4 may be 
separated into two curves. The lower curve sloping upwards 
to the left represents the forenoon data with collector 
efficiency increasing due to a continuous increase in 
radiation during the forenoon. The upper curve sloping 
downward to the right represents the afternoon data with the 
collector efficiency decreasing due to a continuous decrease 
in radiation in the afternoon. The two curves meet at about 
solar noon. In theory, these two curves should be the same 


for a cloudless day. The reasons for the difference in these 
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two curves may be due to: 

1. Energy storage 
The absorbed energy during the forenoon is partly used 
to heat the components of the collector, so the useful 
energy collected is lower and the efficiency is lower . 
In the afternoon, as the solar radiation and ambient 
temperature decrease, some of the stored energy appears 
as useful energy and increases the collector efficiency. 
FiguresoS¢s andaS.Sualsetindzrecatelthetcolhecton 
efficiency continues to increase for 2 to 3 hours after 
solar noon. 

2nveRadrvat ron 
The rate of change of radiation during the forenoon and 
the afternoon may not be the same. Figures 5.3 and 5.5 
indicate that the rate of change in solar radiation 
during the forenoon and the afternoon is not the same. 
This should directly affect the efficiency and the 


(T£i-Ta)/Ht terms. 


Thus, the difference in collector efficiency between 
the forenoon and afternoon depends on: 
a. the mass and heat capacity of the collector 
components 
b. the rate of change in solar radiation. 
The unequal length of the forenoon and afternoon curves 
appears to be related to the (Tfe-Tfi)/Ht term. Some of the 


useful energy used to heat the cold collectors in the early 
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morning is released during the afternoon. This delays the 
lowering of the collector exit air temperature and increases 
the collector operating time. The erratic behavior in 
efficiency at the end of each curve is due to the on-off 


cycling of the system at start-up and shut-down. 


SJo.3 Hourly Efficiency 

in Prqures Sz 1601.5.5, the-plots.of collector 
efficiency over a day suggest that collector approachs 
Steady state at about solar noon because the rate of change 
of storage energy in collector approachs zero. In order to 
investigate this possibility, the records that indicated 
Gontinuousesolar ‘coltection= form the hour thateincluded solar 
noon were selected for study. Ideally, only those days with 
clear sky conditions should be used, but these were few in 
number, So a combination of clear days and partly cloudy 
days were used. A plot of hourly collector efficiency at 
solar noon versus (Tfi-Ta)/Ht is presented in Figure 5.6 and 
the method of least squares was applied to linearize the 
data. The best fit straight line indicates Fr(7ra)e value at 
about 0.54 and the slope indicates the collector heat loss 
coefficient at about 1.58 W/m?-°C. The method was repeated 
using the hour containing solar noon and adding data for the 
previous and following hour. This gave an Fr(ra)e of 0.53 
anc anu ou cr. I 4saw/me— Oo in Figure 5.7 .bTherdata 
exhibits considerable scatter in these two plots, 


particularly for the larger values of (Tfi-Ta)/Ht. The major 
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cause of data scatter is due to clouds that suddenly reduce 
Solar radiation for short periods of time. The (Tfe-Tfi)/Ht 
or (T£i-Ta)/Ht term will increase if solar radiation Ht 
suddenly reduces. But only a small change in (Tfe-Tfi) and 
(T£i-Ta) may occur due to release of energy stored in the 
collectors Therefore, the collector efficiency and 
(T£i-Ta)/Ht terms produce erroneous results. On the other 
hand, the same argument applies for a sudden increase in 
solar radiation as clouds pass. Therefore, an increase or 
decrease in radiation shifts the efficiency plot in opposite 
directions. In addition, the slope represents the heat loss 
coefficient FrU! for the collector and is directly affected 


by data scatter. 


5.5.4 Daily Efficiency 

Agdagpis the=period ofvaycycle! for ‘the collectors” to 
heat up and cool down. Many designers suggest that system 
performance is related to the performance of the collector 
for each day. Thus, the collector efficiency on a daily 
basis was studied. The sum of the useful energy collected, 
the solar radiation and the temperature difference between 
the collector inlet and the ambient temperature were 
calculated on a daily basis. The resulting plot is presented 
in Figure 5.6. unethis plot, the daily ‘data still “shows 
considerable scatter. After the data is linearized by the 
least squares method, the best fit straight line indicates 


an Fr(ra)e of 0.55 and Fru” of 2.56 W/m?-°C. The scatter of 
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daily data may be caused by two reasons: 

1. There may be cloudy periods during the day that reduce 
the solar radiation 

2. There may be a few extremely cold and windy days. This 
results in larger collector losses that lower the 


collector daily efficiency 


5.5.5 Monthly Efficiency 

The plots of instantaneous, hourly and daily basis 
suggest that evaluation over a longer time period may be of 
assistance in reducing the data scatter. A time period of 
one month was selected for study as this period is often 
used in calculations. The solar radiation on the collector 
tilted surface, the useful energy collected, the temperature 
difference between collector inlet and ambient temperature 
were calculated for each month. The resulting plot is 
presented in Figure 5.9 and gives an Fr(ra)e of 0.56 and 
Frul of 2961 Wem? ?C. 

Table 5.8 presents the monthly meteorological data and 
collector thermal efficiency calculated from solar radiation 
and collector air temperature measurements. The normal mean 
monthly ambient temperature and the normal hours of bright 
sunshine are the long-term average monthly values obtained 
from the meteorological records. The R factor is the ratio 
of monthly total radiation on a collector tilted surface to 


monthly total radiation on a horizontal surface. 
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5.6 Summary For Efficiency Methods 
The collector used for module 6 was assembled using 6 

Solaron series 2000 flat plate air collectors by connecting 

two in series and mounting three of these units in parallel. 

This arrangement provided an air heating path length of 

about 3.5 m. The dynamic behavior of the collector makes 

comparison of the actual collector performance to the 
published test performance difficult. Table 5.9 lists the 
collector efficiency calculated for the different time 
periods. The results indicate that all methods produced 

Fr(ra)e values larger than the performance test value of 

0.5. Reasons for this may be: 

1. The performance test value is determined under steady 
state operation with little or no change in energy 
storage in the collector and with little change in solar 
radiation. Under dynamic conditions, energy storage in 
the collector during the morning results ina steep 
Slope for the efficiency curve. During the early 
afternoon, the continuing increase in ambient 
temperature and release of stored energy from the 
cOllectormresubes invian efficiency “curve, that first 
increases and then has a low slope for the rest of the 
afternoon. Representing the data by a best straight line 
results in a larger value for Fr(ra)e. 

2. The accuracy ofthe data may have been, a factor. The 
6rror in mMeasvuming collector flow rate waS about 5 %. 


The error in meaSuring air temperature difference across 
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Table 5.9 Comparison Of Fr(ra)e And Frul 
For Different Methods 


Fr(ra)e Frul 
w/m?2-°C 
Test 0.5 29 
Instantaneous _ —_ 
0.54 La5som # 
Hourly 
Ov53 1r43eAe5 
Daily 0.55 256 
Monthly 0.56 2264 


* Data for the hour including solar noon 
ee Datas for stheshounr=includingssobarsnoonjwithsone hour 
before and one hour thereafter. 
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the collector and air temperature difference between 
room and ambient was about 1.5 %. The solar radiation on 
the collector tilted surface had an error of about 3 %. 
Using these error values, the calculated errors for the 
collector efficiency and the (Tf£i-Ta)/Ht term were 
approximately 6 % and 3 % respectively. This produced a 
7 to 8 % error in the Fr(ra)e parameter. 

The collector heat loss coefficient FrU!t of 1.43 to 
2061°W/m2-°C is lLowerethan®s2.9°W/m?-°C obtained from*the 
rating test. Since insulation had been added to the back of 
the collector, lower values were expected. In addition, the 
Scatter of the data directly affects the FrU+ value. For the 
instantaneous method, energy storage increases during the 
forenoon and decreases in the afternoon. For clear sky 
conditions, the solar radiation follows the same pattern, 
but, when clouds are present, large changes in radiation 
occur, and this can produce ficticious instantaneous 
efficiencies that often very low or very large. For the 
hourly and daily plots, the data scatter indicates the cloud 
effect is present. The monthly FrU! value is probably the 


more representative of the collector heat losses. 
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5.7 Dynamic Behavior Of The Rock Bed Storage Bin 

Although a number of investigators and storage bin 
designers have tried many methods of studying the dynamic 
behaviour of a storage bin, no unique method has been 
proposed. A bin with 2 m* of washed gravel is used for 
energy storage for this thermal system. Due to a shortage of 
data channels, only the air temperature at the top and 
bottom of the rock bed were recorded. This selection allowed 
the calculation of energy added or removed from storage bed. 
In order to study the dynamic behavior of the storage bin, 
plots of air temperature at the top of the storage bin and 
alr temperature at the bottom of storage bin versus the day 
of the week are presented. 

Figures 5.10 through 5.16 show the variation in air 
temperature at the top and bottom of the storage bin during 
a weekly period for each month. The number on the horizontal 
axis indicates the end of the day. Figure 5.10 illustrates a 
period of several days when the bin was overcharged. Figures 
5.11 and 5.13 show the bin behavior during a period of low 
energy storage. Figure 5.14 shows a week when no storage 
occurred. Figures o.12, 5. 1> and i5: 16 ‘show the behavior sof 
the bin when the rock bed temperature wave reached the 
bottom of the bin during a sunny period. 

In Figure 5.15, the third week of February 1982, the 
results for February 15th indicate that there 1s no energy 
storage inside the bin, therefore, the temperature gradient 


across the bin is zero. The plot also indicates that there 
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Figure 5.10 Air Temperature At Top And Bottom Of The 
Storage Bin For The Week Of 2 To 8, 
September 1981 
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Figure 5.10 


a -AIR TEMP. AT TOP GF STORAGE BIN 
+ AIR TEMP. AT BOTTOM GF STORAGE BIN 


DAY 


Alr Temperature At Top And Bottom Of The 
Storage Bin For The Week Of 19 To 25, 
October 1981 
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Figure 5.12 Air Temperature At Top And Bottom Of The 
Storage Bin For The Week Of 5 To 11, 
November 1981 
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Air Temperature At Top And Bottom Of The 
Storage Bin For The Week Of 1 To 7, 
December 1981 
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Figure 5.14 
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Air Temperature At Top And Bottom Of The 
Storage Bin For The Week Of 1 To 7, 
January 1982 
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Figure 5.15 


16 17 18 19 20 


Air Temperature At Top And Bottom Of The 
Storage Bin For The Week Of 15 To 21, 
February 1982 
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-RAIR TEMP. AT TOP GF STORAGE BIN 
-AIR TEMP. AT BOTTOM GF STORAGE BIN 
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Air Temperature At Top And Bottom Of The 
Storage Bin For The Week Of 23 To 29, 
March 1982 
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is little solar energy stored in the bin during days 16 and 
18 of February. This may be due to cloud effect or that any 
solar energy collected was used immediately for space 
heating. Figure 5.15 shows that there are two peaks for days 
17 and 19 through 21 of February. In order to have a more 
detailed study of dynamic behaviour of the rock bin, a plot 
of air temperature at the top and bottom of the bin, 
collector inlet and outlet temperature, and ambient 
temperature versus the time of day for 2 consecutive days, 
20 and 21 of February of 1982, is presented in Figure 5.17 
The system operation modes are also shown on the top of the 
plot. Figures 5.15 and 5.17 indicate that the maximum air 
temperature at the top of the bin occurred shortly after 
noon, decreased to about 40 °C, increased again to about 

48 °C when the system switched to the heat-from-storage mode 
in the late afternoon, and then continually decreased until 
solar collection began the next morning. The air temperature 
at the bottom of the bin remains nearly constant at about 
room temperature except a few degrees increase in the 
afternoon. 

During the forenoon, the temperature of the air 
entering the bin continuouly increased and a temperature 
gradient was established throughout the bin. The top of the 
rock bin would be approaching the air temperature. Shortly 
after solar noon, the air delivered to the bin was at a 
lower temperature than the top of the rock bin. This air 


would start to cool the top of the rock bed, and while no 
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energy was lost, the energy was carried to lower levels of 
the bin. In this way, a temperature wavefront propagated 
Slowly through the bin. Eventually, this wave would reach 
the bottom of the bin and increase the air temperature at 
the bottom of the bin. 

Temperature data presented in Table A-1 (Appendix A-1) 
shows that space heating was required at about 1700 hour. 
The result shows that no energy was available from the 
collectors but energy was available in the storage bin, 
therefore the system operation mode switched from the heat 
storage mode to the heat-from-storage mode with air flow in 
the upward direction. The air temperature at the top of the 
bin increased due to the high temperature in the centre 
portion of the bin and decreased as the stored energy was 
used. Therefore, the first peak that appeared shortly after 
solar noon was due to solar energy sent directly to storage, 
and the second peak that appeared in the evening was due to 
the storage energy in the bin during space heating. 

In Figure 5.11, the first week of September 1981, this 
plot shows that the bin was very hot during days 5 through 8 
of September. Also, the peaks appearing in each of these 
days are different» from days 20:and 21 of February 1982. 
During days 5 through 8 of September, the first peak 
appearing in the early morning was due to the storage energy 
in the bin being used for during space heating , and the 
second peak appearing at about solar noon was due to solar 


energy during energy storage. 
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5.8 Performance Of Storage Bin 

The storage bin is one of the principal components of 
the thermal air heating system, therefore the performance of 
storage bin should directly affect the performance of the 
thermal system. The heat transfers from the walls of the bin 
should be minimized not only to increase the performance of 
the thermal system but also to avoid overheating the module 
in the summer. The rock bed temperature stratification and 
behavior of the bin components could not be studied with the 
available instrumentation. Instead, the monthly energy flows 
to and from the bin were calculated to estimate the energy 
ess for ithe thin. 

Table 5.10 shows that the monthly energy transfers and 
the losses from the storage bin through air leakage or wall 
conduction. During the heating season from September 1981 to 
March 1982, 834.6 MJ of energy were lost from the storage 
bin through wall-conductiom or air leakage. This was) 23.6 % 
of the total solar energy stored in the bin. The losses are 
high but still acceptable. It is interesting to note that 
there are little loss in the months of November 1981 through 
January 1982. This is due to little solar energy stored in 
the bin and the temperature inside the bin is nearly equal 
to the space(room) temperature. There are 34.8 % and 30.7 % 
energy losses in the months of September 1981 and March 1982 
respectively. 

During the heating season, the energy lost through the 


walls and top of the bin contributed to heating the module. 
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Table 5.10 Monthly Energy Loss From The Storage Bin 


Month Yr 
Sep 8 1 
Oct 8 1 
Nov 8 1 
Dec 8 1 
Jan 82 
Feb 82 
Mar 82 
Rotas 
Season 


Monthly 
Energy 
Storage 


Oscol 
(MJ) 


757.4 
364.4 
546.0 
145.5 

44.8 
At9E5 


1247.7 


BOO 


Monthly 


Heat 
From 
Storage 
QOshfs 
(MJ) 


467.0 


Monthly 
Storage 
With 
Electric 
Qseh 
(MJ) 
ZO.6 
37 0 
285.0 
143.8 
42.9 
ees wey 


SOA 


Monthly Percent 
Storage Loss Of 
Bin Stored 
Loss Energy 
Osloss LS 
(MJ) (%) 
26322 34.8 
L588 Zitac 
207.0 a7 
re 8. eee 
29 4,2 
SSa7 Zea 
388.5 2057 
834.6 
236 


Monthly Qscol = Monthly energy entering the storage bin from 
collector during heat-collection mode 
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The rock bed was constructed on the basement floor and the 
lower air plenum was insulated from the basement floor by 
50.8 mm of rigid insulation. Since the measured average 
hourly air temperature for the lower plenum was about equal 
to room temperature, the loss through the bottom of the bin 
would be equivalent to the loss of basement floor insulated 
with 50.8 mm of rigid insulation. Measured basement floor 
losses were about 4 W/m? and remained almost constant 
throughout the year. Correction for the added insulation 
gave an eStimated energy loss of about 7.5 MJ/month. 

Table 5.11 shows that the monthly energy losses from 
the storage bin including the losses of the duct section 
between the collectors and storage bin. A total of 1805.8 MJ 
which represents about 39.9 % of the total solar energy 
delivered to the bin was lost or transfered to the space 
during this heating season. By comparing with Table 5.10 , 
the results show that the losses or energy transfers in 
September 1981 increased from 34.8 % to 61.7 %. It was not 
possible to separate the indoor and outdoor duct losses, but 
comparison with Table 5.10 indicates that the monthly duct 
losses for October to March varied from about 68. to 22 2% of 
the collected energy that was delivered to storage. The 
larger values occurred during December and January. During 
cold months, an increase inolosstoccurs efornhtheyoutdoor 
ducts, A large loss for September [98leresults due to the 
storage bin being overcharged and high collector outlet 


temperature occurring during the heat collection period. 
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Table 5.11 Monthly Energy Loss From Storage Bin 
Including The Duct Loss 


Month Yr Monthly Monthly Monthly Monthly Percent 
Collected Energy Storage Storage Of 
Energy To From With Loss Storage Solar 
Storage Storage Electric Loss Loss 
Qccol Qshfs Oseh Osloss 
(MJ) (MJ) (MJ) (MJ ) (%) (%) 
Sep 8 1 P2098 467.0 28.6 (258i 64.7 54.0 
Ocik 8 1 415.1 oon. 6 SH tis 2 ZOe0 ees 
Nov 8 1 6G HS 25315 20560 TT 8 1266 Theos, 
Dec 8 1 iene Dad 143.8 43.3 23881 Da 0) 
Jan 82 S20 0-0 42.9 14,1 24.7 Deel 
Feb eZ 490.7 152n0 WS 2 boee 9 32n6 108 
Mar 82 1469.7 B6250 bOvaY 604.5 28.4 28.4 
Total A525.2 1805.8 
Season SiS iS, Pad Ve 


Monthly Qccol = Monthly energy delivered from collectors to 
storage bin during heat-collection mode 


Percent of storage loss to the 
total collected energy delivered 
to the bin 
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Percent of Storage Loss 


Percent of storage loss to the 
total collected energy 
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5.9 Energy Transfers From The Thermal System 

Table 5.12 shows the monthly energy transfer or loss of 
the whole thermal air system. The results indicate that in 
woval 83S. 7MMdphwhiehersyaboute2779 “rok the total solar 
gains, were lost or transferred to the space during the 
heating season from September 1981 to March 1982 . In Table 
5.11, the energy transfers or losses from the bin and the 
duct section between the collector and the bin were about 
2Ai.5 AOL the total solar energy collected. The losses from 
the other components of the thermal air system were 
comparatively small. 

Only a small percentage of solar energy was lost during 
the months of November, December and January. The majority 
of the solar energy collected was immediately used for space 
heating and little solar energy was stored in the storage 
bin. For the other months, the energy loss increased as the 


energy delivered to storage increased. 
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Month Yr 


Sep 8 1 
Oct 8 1 
Nov 8 1 
Dec 8 1 
Jan 82 
Feb 82 
Mar 82 
Poca 


Season 


2 Monthly Energy Transfers From The 
Whole Thermal Air System 
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Monthly Monthly Monthly Monthly Monthly Percent 
Collected Heating 
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6. Conclusions And Recommendations 


6.1 Conclusions 


The performance characteristics and the energy flows 


for the air solar heating system for Module 6 at the Alberta 


Home Heating Research Facility have been measured for the 


September 1981 to March 1982 heating period. Analysis of the 


data gave the following conclusions. 


ie 


Theemonthiy solanleontributioniwasefoundsto be 11.8 % by 
the solar required method, 15.7 % by the solar supply 
method) o13.8e%ebynthe burlding¢theathbhossimethod, (d4dm3 % 
by the power records method, and 17.8 % by the Y-Charet 
method. By comparing the values 11.8 % and 15.7 % of the 
direct method, the most probable value of solar 
contribution is about 14 % during the heating season. 
The results indicate that: 

a. The indirect methods are acceptable methods to 
assess the solar contribution. 

b. Long-term monthly solar assessment by the indirect 
methods should be closer to the actual measured 
solar contribution. 

Fr(tra)e values obtained from the collector efficiency 

eurve were 0056randve!54ét onthe thourtyobasisy20.55 on 

the daily basis and 0.56 on the monthly basis. These are 

higher than the performance test value of 0.5. 

The collector heat loss coefficient FrU! was found to be 


1ea3ew/m?-2Gkand MISStw7misf Ceonrtheshourdyibasis, 
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2.56 W/m’?-°C on the daily basis and 2.61 W/m?-°C on the 
monthly basis. These are lower than 2.9 W/m?-°C obtained 
from the performance test. The monthly Fru! value of 
2.61 W/m?-°C is probably the most representative of the 
collector heat losses. 

The monthly collector operating time was found to range 

frome60e2 hoursyotwdanuary ?982rtosi72¢8chourshof sMarch 

1982. The operation hours for each month were about the 

Same as the monthly bright sunshine hours. 

On a cloudless day, the collector thermal efficiency 

curve for the morning was found to be different from 

that in the afternoon. The main reasons are: 

a. The collected energy absorbed by the collector 
components in the morning and released in the 
afternoon. 

b. The rate of change in solar radiation in the morning 
is different from that in the afternoon. 

On a cloudless day, the air temperature and the rock bed 

temperature at the top of the bin reached a maximum 

shortly after solar noon. As the air temperature at the 
collector outlet decreased in the afternoon, the top of 
the rock bed cooled and the maximum rock bed temperature 
moved to a lower level of the bin. Temperature 
stratification in the bin decreased and lower 
temperatures were available when the system changed to 
the heat-from-storage mode. 


There are 3525.3 MJ of solar energy stored in the bin. A 
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total of 634.6 MJ whiGhiis@about 23.6 % of ‘this stored 
energy, was lost from the bin. This energy loss was 
great during the heating season but still acceptable. 
For the whole thermal system, there was 8415.1 MJ of the 
total solar energy collected by the collectors. 
Therefore) 1838.7 MJ, which is about 21.9% of “the ‘total 
collected energy, was lost or transferred to the space. 
These energy losses contributed to the heating load, 
and, during September 1981, resulted in overheating the 


module. 


Recommendations 

The dynamic behavior of the rock bed storage bin can be 
Studied in more detail by installing thermocouples at 
different levels of the bin. 

A cloudless period in August may be used to evaluate the 
adequacy of the size of storage bin for the summer 
period, and to study the insulation requirements to 
prevent overheating of the space due to energy transfer 
from storage bin and ducts. 

The experimental dark test recommended by Edwards and 
Rhee[43] can be used to correlate the instantaneous 
morning and afternoon data for estimating the two 
parameters of Fr(ra)e and Frul, 

Simulation of solar components by existing computer 
programs may be used to study the theoretical 


hour-by-hour performance of the thermal air system. 
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Appendix A-2 


Table A-2 UA Factors Of Module 5 And Module 6 For Some 


Month Yr 
(ND) 
Sep SA 
(0) 
Oct 8 1 
(34) 
Nov 8 1 
(a7) 
Dec 8 1 
7) 
Jan 82 
(9) 
Feb 82 
oy) 
Mar 82 
ey) 
Average 


Electric Power 


Usage 
P5 P6 
(kW) (kW) 
17S8eal 159.4 
So aoe 484.9 
664.3 652.6 


462.0 436.0 


Temperature 
Difference 
of Room and 


Ambient 
DTRA5 DTRA6 
GG) CoG) 


Days With Little Or No Solar Radiation 


Building Heat 
Loss 
Coefficient 

(UA Factor) 


UAS 
(W/°C) 


Voor 


ND - Number of days selected to be investigated and 


From Equation 23729 


Fua 


with little solar radiation 


Monthly Average UA6 


Monthly Average UA5 


SO 


Sie 7) 


Siegen) 
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UA6 
(W/°C) 
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Appendix A-3 


UA Factor Of Module 6 Using January 1982 
As The Analysis Period 


Monthly Electrical Energy Usage 


In January 1982 = 13264.5 MJ 
Monthly Total Temperature Difference 
Between The Room And Ambient = 526256 «-C 


Temperature In January 1982 


From Equation? 3°18 


13264.5 Sai Sher 
= x =—____—— 
26256'..0 3600 

are] 308 AAW oe 
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Appendix A-4 


Conversion Of The Test Performance 


Gollecror Area’ = "9.500 fo) = 1281 m* 


2 CFM/£t? 
0.0102 ne /s-m? 
0.0184 m*/s 


Air Flow Rate 


The original parameters of the test performance 


Pollbv=' 0267 BIU/it? shack = 3.8 W/me—°Cc 


Be we) so =c0 065 


Cp 10060 /,kge © 


pete 70 


pa= Ae coe kos: = 21 


ae 
aC 


From reference of Beckman and Duffie[26] 


inCp 1.204(0.0184) 1006 


Ac g Vee S 


The Standard parameters become 


mCp 
Ac 
Frut = Fout[ —————— ]] 
mCp 
+ FoUL 
Ac 
2538 
= 3.8 [ —————__ ] 
UPAR Sia a ty Als) 
mCp 
Ac 
ol enh aI week 
mCp 
+ FoUL 
Ac 
12238 
= 0.65 [| —————— ] 
( PARSER er SIRS 
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